Discovery of twin kHz QPOs in the peculiar X-ray binary 

Circinus X-1 
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ABSTRACT 

We report the discovery with the RXTE/PCA of twin kHz QPOs in the pccuHar X-ray binary 
Circinus X-1. At eleven different epochs we observed two simultaneous kHz QPOs with centroid 
frequencies of up to ~225 and ^500 Hz and significances of up to 6.3 and 5.5a- for the lower and 
upper kHz QPO, respectively. The frequency range over which the twin kHz QPOs are seen is 
for the most part well below that of other sources, and extends down to centroid frequencies of 
~56 Hz and ~230 Hz respectively for the lower and simultaneously observed upper kHz QPO. 
The strongly variable QPO frequencies and their observed correlations clearly indicate that the 
twin peaks are the kHz QPOs such as typically seen from low-magnetic field neutron stars, and 
not black-hole high-frequency QPOs. This confirms that Cir X-1 is a neutron star, as suspected 
since the detection of Type-I X-ray bursts from the field of the source 20 years ago. The kHz QPO 
peak separation varies over the largest range yet seen, ~175-340 Hz, and increases as a function 
of kHz QPO frequency. This is contrary to what has been observed in other sources but in good 
agreement with predictions of the relativistic precession model and Alfven wave models at these 
low QPO frequencies. Current beat-frequency models would require further modification in order 
to accommodate this. A total of 67 observations showed only a lower kHz QPO; this QPO can 
be followed down all the way to a centroid frequency of ~12 Hz (Q ^^0.5). Two observations 
showed only an upper kHz QPO near 450 Hz. In addition, a strong low-frequency quasi-periodic 
oscillation (LP QPO) is seen in the 1-30 Hz range, as well as two components above and below 
this QPO. The frequency- frequency correlations between the kHz QPOs and the LP QPO are 
in good agreement with those found previously in Z sources, confirming that Cir X-1 may be a 
peculiar Z source, and extend these correlations to frequencies a factor 1.5-2.3 lower. We suggest 
that the low frequency range over which the kHz QPOs occur in Cir X-1 and to a lesser extent in 
(other) Z sources, might be due to a stronger radial accretion flow relative to the disk flow than 
in other kHz QPO sources, possibly related to the nature of the companion star. 

Subject headings: X-rays: binaries. X-rays: individual (Circinus X-1) 
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1. Introduction 

Low-mass X-ray binaries (LMXBs) often ex- 
hibit quasi-periodic oscillations (QPOs) in their 
X-ray flux (van der Klis 2006 for a review). The 
properties of these QPOs vary as a function of the 
location of the source in the track it traces out in 
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X-ray color-color and hardness-intensity diagrams 
(CD and HIDs). Among the low magnetic-field 
neutron-star systems, systematic differences in the 
CD/HID tracks and in the associated rapid X-ray 
variability allow to distinguish two subclasses: Z 
and atoll sources, named after the shapes of their 
tracks (Hasinger & van dcr Klis 1989). Both in 
Z sources and in atoll sources, the QPO frequency 
tends to increase in the same direction as the mass 
accretion rate is inferred to be increasing, at least 
on short time scales (see van der Klis 2001). 

In about 20 LMXBs containing a low-magnetic 
field neutron star, two kHz QPOs have been seen 
simultaneously (twin kHz QPOs) with a centroid 
frequency separation Az/ of a few hundred Hz that 
decreases typically by a few tens of Hz when the 
QPO frequencies themselves increase by hundreds 
of Hz (roughly, from 200 to 900 Hz and from 500 
to 1200 Hz, respectively). If the spin frequency 
of the neutron star, Ug, is known, then Av has 
always been found to be either roughly equal to 
fg, or to J^s/2, depending on the source. The fre- 
quencies of the twin kHz QPOs (called the lower 
and the upper kHz QPO, or Lg and Lu, in order 
of increasing frequency) are tightly correlated to 
each other, following a relation that is the same 
across all sources (Psaltis et al. 1998; Psaltis, Bel- 
loni & van der Klis 1999; van Straaten ct al. 2005). 
Both kHz QPO frequencies also (e.g., Psaltis, Bel- 
loni & van der Klis 1999) correlate with the fre- 
quencies of other variability components such as, 
in atoll sources, the low- frequency quasi-periodic 
oscillation (LF QPO; a sharp QPO peak in the 
1-60 Hz range designated L^p - cf. van Straaten 
et al. 2005 for some of the terminology) and the 
break component (a broad noise component with 
a frequency in the 0.01 20 Hz range designated 
Lb). The components most similar to these in 
the Z sources (and suspected to be physically the 
same) are the horizontal branch oscillation (HBO) 
and the low frequency noise (LFN), but the cor- 
relations of the frequencies of these components 
with those of the kHz QPOs differ somewhat (van 
Straaten ct al. 2003). 

Twin kHz QPOs are a neutron-star signature. 
In some LMXBs containing a black hole, high- 
frequency quasi-periodic oscillations (HP QPOs) 
have been observed. These are different from kHz 
QPOs: they tend to recur near fixed frequencies 
that are different in each source, and when more 



than one HP QPO occurs in the same source the 
QPOs usually appear at different times and at fre- 
quencies often roughly consistent with a fixed 2:3 
ratio (Strohmayer 2001; McClintock & Remillard 
2006). 

Numerous models have been proposed to ex- 
plain the nature of the kHz QPOs and describe 
their correlations. Some refer to a blob orbit- 
ing inside the accretion disk, either self-luminous 
or scattering the X-rays, or producing a hot spot 
on the neutron-star surface by the infall of mat- 
ter from the blob while it is orbiting at the inner 
edge of the disk. These models produce a rotating 
pattern of radiation causing a so-called 'beaming' 
modulation of the observed flux. The orbital fre- 
quency of a particle in an equatorial circular orbit 
around a point mass with mass M and angular 
momentum J (approximating the spacetime met- 
ric around the neutron star by Kerr spacetime) 
is = v^GM/ (47r2r3)/(l + jirg/rf/^), where 
rg = GM/c^ and j = Jc/GM'^. The correspond- 
ing radial cpicyclic frequency in a slightly eccentric 
orbit and for negligible frame dragging (j « 0) is 
i^r ~ z^0a/1 - Qrg/r. 

In the relativistic precession model (Stella & 
Vietri 1999) as well as in the sonic-point spin- 
resonance model (Lamb & Miller 2003) the upper 
kHz QPO Lu is identified with i'^. The lower kHz 
QPO Li is interpreted in the first case as due to 
the periastron precession frequency v^ — Vr, and in 
the second case as a beat between at the sonic 
radius and the spin frequency (at approximately 
Vff, - Vspin or - V spin/ 2 depending on condi- 
tions at the spin-resonance radius where the beat 
emission is generated). Kluzniak & Abramow- 
icz (2001, 2003) suggest that QPOs occurring in 
neutron stars are similar to those seen in black 
holes and result from resonance frequencies in the 
disk, but in neutron stars are tuned away from 
the fixed frequencies seen in black holes by addi- 
tional physical effects (but see also Belloni et al. 
2005). Models related to orbital motion are in 
general more predictive of the frequencies of the 
QPOs than of their amplitudes and coherences, 
as constructing a model for the modulation mech- 
anism or describing the damping of the oscilla- 
tions is more difScult than explaining these fre- 
quencies as a (combination of) general-relativistic 
orbital and epicyclic frequencies. A maximum fre- 
quency according to such models would be the 
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orbital frequency at the innermost (marginally) 
stable circular orbit (ISCO) of the compact ob- 
ject visco » (cV27r63/2GM)(l + 0.75j) for Kerr 
spacctimc or, under certain conditions, the or- 
bital frequency at the marginally bound orbital 
radius. Both these frequencies scale inversely with 
mass. For an extensive review on timing fea- 
tures of LMXBs and QPO models see van der Klis 
(2006). 

Circinus X-1 is a galactic X-ray source located 
in the galactic plane at a distance that has vari- 
ously been reported to lie in the range 4-12 kpc 
(Goss & Mcbold 1977; Glass 1994; Jonker & Nele- 
mans 2004; laria et al. 2005). The system is 
thought to have a 17-d highly eccentric (e ~ 0.8) 
orbit leading to periodically enhanced mass trans- 
fer at periastron (Murdin et al. 1980). From opti- 
cal and infrared observations a sub-giant compan- 
ion star of 3-5 Mq was suggested (Johnston et al. 
1999). More recently, Jonker et al. (2006) found 
evidence for supergiant characteristics in the op- 
tical spectrum of the system which might indicate 
an even higher companion mass. So, Cir X-1 is 
formally a high-mass X-ray binary. However, its 
X-ray phenomenology corresponds to that of a low 
magnetic- field compact object such as more com- 
monly found in LMXBs. The strong radio emis- 
sion (see, e.g., Fender et al. 2005), highly relativis- 
tic jets (Fender et al. 2004), hard X-ray emission 
(laria et al. 2001; Ding et al. 2003) and very strong 
X-ray variability (Jones et al. 1974; Samimi et al. 
1979; Oosterbroek et al. 1995), resemble proper- 
ties generally more typical of black holes. How- 
ever, type-I X-ray bursts seen in 1986 from the 
field of the source and most likely originating in 
Cir X-1 itself (Tennant et al. 1986) classified it 
as a likely neutron star, with later spectroscopic 
analysis (Maisack et al. 1995; Done & Gierlihski 
2003) supporting this at least at a phenomeno- 
logical level. No bursts were reported from the 
source since 1986. Low- frequency noise, a clear 
low- frequency QPO L^p and a broad feature at 
100-200 Hz were seen with EXOSAT (Tennant 
1987). At low flux levels Cir X-1 exhibited atoll- 
source behavior (Oosterbroek et al. 1995) includ- 
ing strong band-limited noise, while subsequent 
observations at higher flux on occasion, and in 
addition to a wide variety of non-Z-like patterns, 
showed a Z-source-like track in the HID (Shirey et 
al. 1999). On the horizontal branch of this track, 



Shirey et al. (1996, 1998) found the LF QPO to 
move between 1.3 Hz and 32 Hz, in strong correla- 
tion with the break frequency at lower frequencies 
and the frequency of the broad high-frequency fea- 
ture, now seen to range from 20-100 Hz. Psaltis, 
Belloni & van der Klis (1999) interpreted this 
broad peak low-frequency, low-coherence ver- 
sion of the lower kHz QPO as seen in Z and 
atoll sources. They showed that its frequency fits 
on the relation seen in those sources between the 
frequencies of and Llf, extending it down to an 
Lfi centroid frequency of ~20 Hz and, interestingly, 
providing a link with black-hole phenomenology. 
No twin kHz QPOs were reported from Cir X-1 
up to now. 

In this paper, we concentrate on the high- 
frequency part of the power spectra of Cir X-1 
and perform an extensive survey for kHz QPOs. 
We find that twin kHz QPOs do in fact occur in 
Cir X-1. We describe our methodology in Sec- 
tion 2, present our results in Section 3 and com- 
ment on consequences for phenomenological sys- 
tematics and theoretical models in Section 4. 

2. Observations and Analysis 

We used the data on Cir X-1 obtained with the 

proportional counter array (PCA; Jahoda et al. 
2006) on board the Rossi X-ray Timing Explorer 
(RXTE) from March 1996 tiU January 2005 hsted 
in Table 1. This includes all data publicly avail- 
able in February 2005 as well as the then propri- 
etary data set 90025, and comprises about 2 Ms 
of data distributed over 21 RXTE programs in 9 
RXTE Observing Cycles: a total of 497 RXTE ob- 
servations each with typically 3-6 ks of useful data 
and identified by a unique RXTE Observation ID. 

The source was observed with 13 different com- 
binations of instrumental modes covering the full 
energy range (channels 249, effective energy 
range 3-60 keV), and one more starting from chan- 
nel 8 (3.7 keV) with which 21 observations were 
taken in 2003 (Cycle 8), when the intensity was 
low (below 0.4 Crab) and from which we eventu- 
ally did not report any power spectra. We divided 
each observation into continuous time segments of 
16, 64, 128 or 256 s, depending on the time res- 
olution such that they contained 2^*^ points. The 
count rates from all of the PCA's available Propor- 
tional Counter Units (PCUs) were summed and 
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Fourier-transformed. The resulting power spectra 
were averaged to give one average power spectrum 
per observation. 

Most of our power spectra (378 of the 497) ei- 
ther showed no appreciable broad-band power at 
all or only a broad plateau-like feature at low fre- 
quency (< 10 Hz). These power spectra were 
inspected visually for high-frequency features, if 
necessary followed by fitting using the procedures 
described below in order to determine the signif- 
icance of any features. Taking into account the 
number of trials inherent in this search, no sig- 
nificant high-frequency features were found. The 
most significant feature (at ^721 Hz and with 
Q^8.5) had a single-trial significance of 3.8it. 
More complex power spectra, containing a signifi- 
cant and narrow {Q >1) LF QPO in the 1-50 Hz 
range and at least one more feature described well 
by a Lorentzian, were seen in the remaining ob- 
servations. We performed a full analysis of these 
119 power spectra using the multi-Lorentzian fit 
method described below. 

Since this work concentrates on kHz QPOs, 
careful subtraction of the deadtime-modified 
Poisson-noise power is required. A full descrip- 
tion of our approach is provided in Appendix A. 
We used the expression proposed by Zhang et al. 
(1995, 1996) to describe the Poisson spectrum, 
but with values of the event deadtime and VLE 
window of td = 8.87 /is and tvLE = 162/is based 
on direct fits of that model to our data. As ex- 
plained in the Appendix, this method is slightly 
more conservative with respect to the significances 
of the twin kHz QPOs than the other methods 
used in previous work, and provides a consistent 
estimate of the deadtime-modified Poisson-noise 
power spectrum throughout our data. 

After subtracting the noise, we renormalized 
the power spectra to source squared fractional rms 
amplitude and fitted them using standard chi- 
squared methods. We used a multi-Lorentzian 
fit function in the {vmax, Q) representation (Bel- 
loni et al. 2002), where the so-called characteris- 
tic frequency Vmax of the Lorentzian is given by 
Vynax = v^i'o + with vq the Lorentzian's cen- 
troid frequency and A its HWHM, and where its 
quality factor is Q = i^o/2A. A power law was 
added at low frequencies when needed. The fits 
were made without any reference to expected fea- 
tures, and include all significant (> Sct; all signif- 



icances quoted herein are single-trial) Lorentzians 
that could be fitted. All errors quoted are 1-cr 
single-parameter (Ax^ = 1). 

To recover the Lorentzian centroid frequency 
vq, we used the expression = VmaxI + 1/4(5^- 
Since all Lorentzians for which we applied this 
conversion are rather sharp-peaked (Q > 1), the 
errors computed for Vjnax are approximately the 
same as those for z/q. For the twin kHz QPOs, 
we explicitly re-estimated the errors by refitting 
in centroid-frequency representation, and typi- 
cally found the errors to differ by <5% (usually 
<1 Hz), and negligible differences in the parame- 
ters themselves. For the other fits we adopted the 
Vraax errors also for the centroid frequencies. In 
the following, all frequencies quoted are in terms 
of the characteristic frequency v^^ax unless oth- 
erwise noted. Subscripts denote the variability 
component {vu for upper kHz QPO, etc.); cen- 
troid frequencies are denoted with superscript 
as follows: v\. 

We used the 16-s time-resolution Standard-2 
mode to calculate X-ray colors as in Altamirano 
et al. (2005) but with bands as defined by Shirey 
et al. (1999). Dropouts were removed, back- 
ground was subtracted and deadtime corrections 
were made. All active PCUs were used for this 
analysis. The energy bands chosen for colors and 
intensity are given in the captions to Figs. 3 and 
4 - energy channels were interpolated to approx- 
imate the exact same bands in all data sets. In 
order to correct for the gain changes as well as the 
differences in effective area between the PCUs, we 
normalized all color and intensity estimates to the 
corresponding Crab values (Kmdkers at al. 1994; 
van Straaten et al. 2003) that are close in time but 
in the same RXTE gain epoch. 

3. Results 

In 80 of the 119 fitted power spectra we found 
a QPO which we interpret as a kHz QPO (see 
below). In 7 of these we found a significant (> 
3ct) second kHz QPO. A second kHz QPO with 
a significance between 2.5 and 3cr was found in 4 
more cases. We discuss whether these are real or 
not in Section 3.1. 

The best power spectrum exhibiting twin kHz 
QPOs (data set H; observation 20094-01-01-01) is 
shown in Fig. 1. Twin kHz QPOs, 6.3 and 5.5cr 
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significant, arc clearly seen near 140 and 460 Hz, 
respectively, as well as a strong LF QPO near 
10 Hz, broad noise, and a power law at low fre- 
quency. The power spectra of all 11 observations 
with twin kHz QPOs are shown in Fig. 2 together 
with all fitted components. The prominent narrow 
feature peaking somewhere between 1 and 30 Hz 
is the LF QPO Llf- It is seen in all 119 observa- 
tions, sometimes as a double peak, which may be 
due to averaging over a moving feature. When due 
to this it was necessary to fit the LF QPO with 
two Lorentzians, we adopted the frequency of the 
most significant (usually, the lower-frequency) of 
the two as that of the LF QPO in further anal- 
ysis. A low- frequency noise or break component 
was seen below i/lf in all but three cases while an 
additional component was detected at or above 
ulf in more than half (78) of the power spectra. 
The nature of this component is further discussed 
in Section 3.3. Power laws were needed in only 
nine cases, and always only produced appreciable 
power below ~1 Hz. 

The average X-ray colors of each observation 
with kHz QPOs are plotted on top of a 16-s time- 
resolution color-color diagram (CD) of all Cir X-1 
data in Fig. 3. The observations with twin kHz 
QPOs and single kHz QPOs are set out by differ- 
ent symbols, as are the 16-s data already reported 
by Shirey et al. (1999). Twin and single kHz QPOs 
occur over a range of colors; the frequency of the 
lower kHz QPO tends to increase towards both 
lower and higher colors, but the correlations are 
weak. A hardness-intensity diagram following the 
same conventions is provided in Fig. 4. Most of the 
observations showing twin kHz QPOs were within 
a quite narrow band of intensities between 2545 
and 2770 c/s/PCU, or about 1.0-1.1 Crab (ex- 
cept data set H at 1991 c/s/PCU or ~0.8 Crab, 
which is also the overall hardest kHz QPO obser- 
vation; the uppermost black triangle in Figs. 3 and 
4). Presumably related to this, we find no clear 
correlation between kHz QPO frequency and in- 
tensity, although there is some weak c;vidc;nce for 
parallel tracks in frequency vs. intensity, as has 
been seen in other sources. Most kHz QPOs, and 
all twins, are found well above the Z tracks iden- 
tified by Shirey et al. (1999) which are confined 
to Crab-normalized broad color < 0.6, i.e., most 
QPOs are associated with a considerably harder 
overall spectrum than that in the Z tracks. Fig. 5 



shows a long term light curve of the source where 
each data point is the average of one observation. 
It can be clearly seen that most of the observations 
with kHz QPOs, and all with twin kHz QPOs, oc- 
curred at relatively early epochs, when the average 
source flux was higher than later in time, and that 
most cluster near intensities of 1 Crab. 

To investigate the relations between the fre- 
quencies of the power-spectral components we plot 
the characteristic frequency of all features vs. ly^p 
in Fig. 6. The points can be easily separated into 
five groups of correlated frequencies. The LF QPO 
{small filled circles) is the most significant feature 
in our power spectra. A dispersed group of points 
below the LF QPO (triangles) corresponds to the 
low-frequency noise or break component (see also 
Section 3.3). At a frequency between 1 and 3 times 
that of the LF QPO we see a group of points as- 
sociated with a component (open circles) that is 
usually dominated by a harmonic to the LF QPO 
(see Section 3.3). Parallel to that and at somewhat 
higher frequencies there is an almost continuous 
series of points extending from ^18 up to ~380 Hz 
(large filed circles). Yet another group of points 
above these has similar characteristics, although 
less dense, and extends from ~235 up to ~580 Hz 
(squares). Prom the high frequencies to which they 
extend, and their frequency variability over more 
than an order of magnitude, we identify these lat- 
ter two components as kHz QPOs. We conclude 
that the highest-frequency group (squares) corre- 
sponds to the upper kHz QPO (L„) and the group 
below that (large filled circles) to the lower kHz 
QPO (L^). These components coincide with 
the lower kHz QPOs as tentatively identified by 
Psaltis, Belloni & van der Klis (1999, see Sec- 
tion 1), and are found here over a wider frequency 
range. Note that in accordance with previous work 
we use the term kHz QPOs for components that 
can have high frequencies, but are also seen down 
to a characteristic frequency as low as (in our case) 
18 Hz. 

3.1. Tv^rin kHz QPOs 

As noted in Section 3, Fig. 2 shows the power 
spectra of the observations exhibiting twin kHz 
QPOs that were both significant at >3(t (up to 
17cr for L( and up to 5.5cr for Lu), as well as those 
with a significance >3cr for one and >2.5cr for the 
other. Prom Figs. 3 and 4 we see that the points of 
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the 2.5-3(7 group fall right among the >3(t group. 
The QPO frequencies of the two groups are also 
in the same range, and the frequency correlations 
seem to hold for both groups of points (Fig. 6). 
We conclude that it is safe to include the four ob- 
servations with a kHz QPO between 2.5 and Str in 
our sample of twin kHz QPOs; of course they have 
somewhat larger statistical error bars. Weak fea- 
tures above i/„ visible in Fig. 2 (e.g. in power spec- 
trum E), and in between Vu and vi> (e.g. in power 
spectrum G) are not significant. We remain thus 
with 11 twin kHz QPOs. Table 2 lists their best- 
fit ccntroid frequencies, coherences and fractional 
rms amplitudes, as well as the peak centroid sepa- 
ration Ai/ and the centroid frequency ratio. Note 
that Ai' tends to increase with QPO frequency 
and that the frequency ratio is inconsistent with 
2:3. This is further discussed in Section 4. 

3.2. Single kHz QPOs 

Apart from the 11 twin kHz QPOs we measured 
an additional 69 single ones. Their frequency 
shows a strong correlation with that of Llf- As 
can be seen in Fig. 6, these single kHz QPOs are 
nearly all consistent with being lower kHz QPOs: 
the 11 lower kHz QPOs among the twin peaks are 
clearly part of the same group of points. Only 
the two highest-frequency single kHz QPO points 
(20094-01-01-00 and 20094-01-01-000) seem to di- 
verge from this correlation (particularly in terms 
of centroid frequencies, not shown). Since these 
points seem instead to extend the relation for the 
upper kHz QPOs, and because a less significant 
feature is present in those power spectra at fre- 
quencies consistent with a lower kHz QPO, we 
classify these two QPOs as L„, with remain- 
ing formally undetected due to statistics. This is 
consistent with the decrease in lower kHz QPO 
amplitude with increasing frequency (Fig. 7). The 
best-fit power law for the lower kHz QPO points 
is vt = (19.3 ±0.6) X Hz. with a x^/dof 

of 17.8 for 76 degrees of freedom, clearly indicat- 
ing the existence of non-statistical scatter in these 
frequencies; here and below the error bars quoted 
for power-law fit parameters were multiplied by 
y/x^/dof to take account of this scatter. For L„, 
we similarly find Vu = (130 ± 30) x zy0|0±0-09 Hz; 

/dof = 3.4 for 11 degrees of freedom. 

The quality factors of the kHz QPOs ranged 
up to Qi^Q for and up to (3~2.5 for L„ but 



were more often around 1. No clear trend was 
detected in the Q values as a function of QPO fre- 
quency. Averaging over a moving feature may in 
some cases have broadened the peak somewhat, 
but experiments with data subsets show this ef- 
fect to be minor down to accessible time scales 
(^^5000 s). Only in two cases evidence was de- 
tected of a deviation from the standard Lorentzian 
profile strong enough to allow a significant addi- 
tional component to be fitted. We used the results 
of a single Lorentzian fit also in these two cases. 

Fig. 7 shows the fractional rms amplitude of 
each kHz QPO as a function of its own charac- 
teristic frequency. Some general trends are seen, 
where the rms of monotonically decreases (but 
more rapidly below a of 200 Hz than above that 
frequency) and the rms of L„ varies much less, but 
may first rise and then fall again. This is some- 
what reminiscent of the amplitude variation of the 
kHz QPOs in the Z source GX 5—1 (Jonker et al. 
2002) and different from what has been observed 
in a number of atoll sources (Mendez et al. 2001). 
Both QPOs appear to vary in amplitude gradually 
and tend to become undetectable when they fall 
below an rms amplitude of ~1.5%. 

3.3. Other power spectral components 

The lower frequency power spectral compo- 
nents are not the prime objective of this work, but 
we briefly report here on a number of immediate 
findings concerning these components that emerge 
from our analysis. 

The component above L^p often is rather co- 
herent (with Q up to ^8), and when it is (for 
Q > 2), it has a ccntroid frequency twice that of 
the LF QPO (Fig. 8) to good precision, indicating 
that mostly this component is the harmonic L2LF 
of the LF QPO. However, when it is less coherent 
the component often appears to be dominated by 
broad structure around or just above the LF QPO. 
It can then have frequencies down to a character- 
istic frequency as low as that of Llf (see Fig. 6), 
and to well below that component in terms of cen- 
troid frequency (see Fig. 8). 

Although an increase in the break frequency 
of the low-frequency noise with v^f is evident in 
Fig. 6, the scatter is large. The power spectra 
at low frequencies may be more complex than al- 
lowed in our fits, which mainly concentrate on kHz 
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QPOs. Possibly, two parallel groups of points arc 
present, indicating two break frequency compo- 
nents as observed in other sources {Lh and Lb2; 
cf. van Straatcn ct al. 2003). This seems for ex- 
ample to be the case in power spectrum F. More 
work is needed to resolve this issue. 

4. Discussion 

We have discovered twin kHz QPOs in Cir X-1 
in 11 separate RXTE/PCA observations. Among 
these 11 observations the centroid frequency of 
the upper kHz QPO ranged from 229±18 to 
505±51 Hz and that of the lower kHz QPO from 
56.1±1.3 to 226±18 Hz. The lowest-frequency 
twin kHz QPOs observed so far, in the accret- 
ing millisecond pulsar XTE J1807-294, had ~127 
and ~353 Hz centroid frequencies (Linares et al. 
2005); our Cir X-1 data extend to well below this. 
A total of 80 observations showed at least one kHz 
QPO, with centroid frequencies extending down 
to as low as 12.5 Hz. The centroid frequency sep- 
aration of the twin kHz QPOs varies by a factor 
of two, spanning a range, 173±18 to 340±64 Hz, 
similar to that covered by all previous twin kHz 
QPO detections in other sources taken together. 

In order to further compare the kHz QPOs in 
Cir X-1 with those in other sources we examine the 
locus of these QPOs in two frequency-frequency 
correlation diagrams that have been considered in 
previous works. In Fig. 9 we plot v'^ and jyip vs. 

in the Psaltis, Belloni & van der Klis (1999) 
diagram containing data for both neutron stars 
and black holes. The lower kHz QPO points of 
Cir X-1 had previously been used by these authors 
to suggest a link between the lower kHz QPO in 
Z and atoll sources above a of 200 Hz, and 
a broad feature in black holes and weak X-ray 
burst sources below a of 10 Hz (both plotted 
in Fig. 9 using open circles). This latter broad 
feature was dubbed L^ow by van Straaten et al. 
(2005), who based on the relations of amplitude 
and coherence with frequency, and on systematic 
shifts in the frequency-frequency correlations ob- 
served in some millisecond pulsars questioned the 
identification with (van Straaten et al. 2003, 
2005). Our identification of the twin kHz QPOs in 
Cir X-1 confirms that at least in Cir X-1 this fea- 
tTire is indeed the lower kHz QPO as proposed by 
Psaltis, Belloni & van der Klis (1999). Our newly 



discovered upper kHz QPO points smoothly ex- 
tend the i/^ vs. relation to frequencies a factor 
4 lower than originally reported (Psaltis, Belloni 
& van der Klis 1999; upper track in Fig. 9). The 
best- fit power law to all points including Cir X-1 is 

= (23±3) xi/0^°-^^^°-°^^ Hz {xVdof = 20 for 82 
degrees of freedom). As can be seen in Fig. 9, this 
relation provides an excellent match to Cir X-1 as 
well as to the data on the other sources. 

In Fig. 10 we plot the characteristic frequencies 
of all other power spectral components in the 11 
twin kHz QPO power spectra vs. in a diagram 
originating with van Straaten et al. (2002) con- 
taining data for atoll sources, weak burst sources 
and Z sources. In this diagram it can be seen 
that the vs. Uu relation, as smoothly extended 
with our Cir X-1 points, does not connect to the 
i^iow vs. Vu relation, supporting van Straaten et 
al.'s suspicion that and Lfow may be different 
components. Further work on the other sources 
contributing to these correlations is required to 
definitively resolve this issue; this is beyond the 
scope of this paper. 

The large range over which the twin kHz QPO 
frequencies vary, the strong deviation from a 2:3 
ratio and the frequent simultaneous appearance of 
the two peaks are very different from what is seen 
in black hole high-frequency QPOs. Instead, the 
way in which the kHz QPO frequencies correlate to 
each other and to the frequencies of other variabil- 
ity components is consistent with what is seen in 
neutron stars. These observations therefore sup- 
port the identification of Cir X-1 as a neutron star 
(cf. Section 1). In Fig. 10 it can be seen that the 
LF QPO in Cir X-1 fits well on the Lh track in 
this diagram, but the low-frequency noise or break 
component does not match well to the Lb track. 
Cir X-1 shares both these characteristics with the 
Z sources (van Straaten et al. 2003), giving further 
support to the proposal of Shirey et al. (1999) that 
Cir X-1 is a (peculiar) Z source. The low coher- 
ence, relative weakness, and low frequency range 
of our twin kHz QPOs as well as the detection of 
a LF QPO harmonic (such as also commonly seen 
in Z source HBOs) are also in accordance with a 
Z source classification. However, the band-limited 
noise can be much stronger (Oosterbroek et al. 
1995), and the kHz QPOs can reach much lower 
frequencies in Cir X-1 than has been observed in 
any (other) Z source. 
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The extended range of kHz QPO frequencies in 
Cir X-1 allows unique tests of theoretical models. 
Perhaps most strikingly, the increase in the fre- 
quency separation [s.v with upper kHz QPO ccn- 
troid frequency v^, which has only been observed 
previously, somewhat marginally, in GX 17+2 
(Homan et al. 2002) and on a single occasion in 
4U 1728-34 (Mighari et al. 2003), is in accor- 
dance with a prediction of the relativistic preces- 
sion model (Stella & Victri 1999) for the behavior 
of Aj/ at low QPO frequency. 

In Fig. 11 we show a plot of the centroid fre- 
quency separation = v'^ — vs. v''^. From 
an F-test for additional term we find that in 
our data, Az/ is inconsistent with being constant 
at the 5.2(7 level; a linear fit yields a line with 
slope 0.6±0.1. The relation predicted by the 
precession model for negligible frame dragging, 

Az/ = z/O-y/l - 6 /c^ describes our 

data points well {drawn curve). The best fit to just 
our Cir X-1 points jdof = 0.2) gives a mass for 
the neutron star of M = (2.2 ± 0.3)Mq. This is 
high, but agrees with inferences drawn from this 
model for other sources (Stella & Vietri 1999), 
and indeed the best-fit curve passes through the 
cloud of scattered points corresponding to all other 
sources (which of course can be expected to have 
different masses). Note, however, that as already 
discussed by Stella & Vietri (1999), for some other 
individual sources no exact match to this model is 
possible without introducing additional physical 
assumptions. We further note, that Alfvcn wave 
oscillation models (Zhang 2004; Rezania & Sam- 
son 2005; Zhang et al. 2005) which also predict Ai/ 
to increase at low t'", can fit our data as well. 

As is evident from the data presented in the 
rightmost column of Table 2, while the twin kHz 
QPO frequencies are clearly inconsistent with a 2:3 
ratio, a 1:3 ratio is closer, although still formally 
inconsistent with the data at 99% confidence. A 
1:3 ratio corresponds to one of the resonances that 
have been considered by Kluzniak & Abramow- 
icz (2003). However, since that ratio would arise 
from a combination frequency of the main 2:3 reso- 
nance, the weakness of the main resonance (which 
is not significantly detected) would in such an in- 
terpretation remain unexplained. Beat-frequency 
models, based on a constant frequency separation 
Az/ related to the spin frequency, would have diffi- 



culties explaining the factor 2 change in frequency 
separation of the kHz QPOs of Cir X-1. The mod- 
ified beat-frequency model proposed by Lamb & 
Miller (2001) to explain the modest decrease of 
Ai/ with QPO frequency previously observed in a 
number of sources would need further modifica- 
tion to explain the considerable increase that we 
observe instead. The range of upper kHz QPO 
frequencies we observe (229-505 Hz) is nearly en- 
tirely below that expected in the sonic-point beat- 
frequency model (500-1200 Hz; Miller et al. 1998). 

Although our discovery of twin kHz QPOs 
in Cir X-1, as well as the observed frequency- 
frequency correlations, put this system on a foot- 
ing with the neutron-star LMXBs, the fact that 
frequency range over which the upper kHz QPOs 
are observed is mostly well below that of other 
sources further confirms the exceptional nature of 
the source. These low frequencies might be taken 
to indicate the presence of a massive neutron star, 
if for example the upper kHz QPO would de- 
rive from the orbital frequency at some radius set 
by general relativity, scaling like 1/M. However, 
when taken at face value the fact that the observed 
frequencies are strongly variable makes this inter- 
pretation unlikely. The possibility remains that 
the highest kHz QPO frequencies in a source de- 
rive through general relativity from M, and the 
observed range of frequencies in turn derives from 
this maximum. 

Most proposed models interpret the strong vari- 
ations of the frequencies as due to variations in 
the inner radius of the Keplerian disk flow. This 
radius could be set by a number of different mech- 
anisms, such as electromagnetic stresses (in which 
case the inner radius would be the radius of the 
magnetospliere), or alternatively radiative stresses 
(radial radiation force and radiation drag on the 
accreting matter; Miller et al. 1998). Within the 
framework of such ideas the low kHz QPO fre- 
quencies in Cir X-1 would then imply a relatively 
stronger magnetic field for the neutron star, or al- 
ternatively a larger ratio of radiation flux density 
to matter density at the inner disk edge. In the 
description proposed by van der KUs (2001) for the 
kHz QPO 'parallel tracks' phenomenon, relatively 
larger radiative stresses might occur if a larger 
fraction of the accretion takes place not through 
the (thin) disk, but in a more radial inflow, lead- 
ing to relatively more radiation, and hence more 
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radiative stress, being produced at the inner disk 
edge per gram of accreting disk material. The ex- 
istence of such a stronger radial inflow might pos- 
sibly be related to the unusual systemic character- 
istics of Cir X-1, such as the periodic episodes of 
rapid mass transfer and/or the large companion 
mass and likely, radius. 

It is interesting to note, that a tendency to- 
wards relatively low kHz QPO frequencies is al- 
ready known to be one of the characteristics dis- 
tinguishing Z sources from other neutron star 
LMXBs. As noted above, Cir X-1 shares several 
more characteristics with the Z sources. So, Cir X- 
1 may possess the same physical property that 
makes Z sources stand out, but more extreme. If 
that property would be a stronger magnetic field, 
as has been discussed before (e.g., Psaltis et al. 
1995), then additional ingredients (such as smear- 
ing of rapid variability by scattering in material 
surrounding the neutron star) are required to ex- 
plain why no coherent pulsations are observed. If 
instead a stronger (radial flow / disk flow) ratio 
causes the Z-source phenomenon (possibly, like in 
Cir X-1, related to systemic characteristics such 
as a larger companion radius), then the radiative 
transfer of the emerging radiation through this in- 
flow may be what causes the spectral effects (the 
Z track), and the larger inner disk radius what 
causes some of the timing effects (lower kHz QPO 
frequencies) that set these sources apart. 
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A. The determination of the Poisson noise spectrum 

For a high-frequency analysis an accurate description of the deadtirnc-modified Poisson noise power spec- 
trum is required. In the RXTE/PCA this spectrum is mainly determined by the values of the event deadtime 
td and the VLE- window^ tvLE- The technical descriptions of RXTE give only approximate values for these 
parameters, td~10 and tvLE=^^0 M^j respectively (Zhang et al. 1996); values of 8.5/is have also been used 
for td (Revnivtsev et al. 2000, Klein- Wolt et al. 2004, Wei 2006), while a method using additional small shifts 
was introduced by Klein- Wolt et al. (2004) to compensate for deviation from the observations. Recently, a 
calibrated VLE window value of 170/is was reported by Jahoda et al. (2006). Since these values could in 
principle vary and in some cases depend on source spectrum, and since the values mentioned above do not 
agree with our observations to high precision (see Fig. 12), we instead estimated the deadtime parameters 
directly, making good use of the large sample of data available for Cir X-1. 

We used 84 representative RXTE observations of Cir X-1 distributed over all epochs and with count 
rates between 300 and 6000 c/s/PCU. No dependencies of the deadtime parameters on epoch or count rate 
were found. Observations with high time resolution were preferred but lower-resolution ones, as well as 
observations with extremely high resolution (Nyquist frequency ~33 kHz) all gave consistent results. In 
order to avoid power contributions from known source variability components, we looked only at frequencies 
above 1.6 kHz (up to the Nyquist frequency). All observations used for this analysis covered the whole 
energy range (channels 249). 

We used the function approximately describing the Poisson power spectrum according to Zhang et al. 
(1995, 1996): 

P. = 2- ,rotd (l - ^) - ^^rotd (I) cos2..i. + 2r...rotk. (^|^) ^ ' 

where td and tvLE were defined above, N is the number of points in the time series, ro is the count rate per 
detector, Vvle is the VLE count rate per detector, and tb is the time bin width, so that l/2tb is the Nyquist 
frequency. By using td and tvLE 8^ free parameters and fixing the other ones to their known values, we fitted 
this function to the power spectra above 1.6 kHz for each of the 84 observations. All fits represented the 
data well, without need for shifting, for a td value spanning only a narrow range, while tvLE was also well 
determined. Based on these measurements of the deadtime parameters and their distributions we adopted 
the values td=8.87/us and tvL£:=162/is to use throughout our entire analysis. These numbers are in the range 
of values given by the RXTE team (Jahoda et al. 2006) and used by various other authors (Revnivtsev et al. 
2000; Wei 2006). Whether the above values accurately describe the deadtime values for other observations 
needs to be tested using sources with other spectral shapes and flux levels, as well as for the other VLE 
settings. 

As a test we compared our fit results to those we would have obtained with other choices for the deadtime 
modified Poisson spectnmi calcidation. Deadtime parameters td=10fis, iv'LE=170/iS gave a very bad match 
to many of our power spectra and in those cases produced no reasonable fits, but some previously used 
deadtime parameter values (near 8.5, 150/is) gave reasonable results. We refitted the power spectra for the 
observations showing twin kHz QPOs (Section 3.1) using 8.5 and 150/ks and found no systematic deviations 
in the kHz QPO frequencies and coherences compared to fitting with our adopted values. This was also the 
case when using the method of shifting the Poisson noise spectrum to better match the power spectra at 
high frequencies. We did find that our adopted method tended to be slightly more conservative with respect 
to detecting high-frequency features than the shifting method, where occasionally an extra component was 
necessary at high frequencies which was not significant with our method (e.g., power spectrum E of Fig. 2). 
In any case, twin kHz QPOs reported as significant herein were also significant with these two previously 
used methods of estimating the deadtime modified Poisson noise spectrum. 



^In our data only the 3rd VLE window setting appears (a.k.a. setting 2; ~150 ^is) and for that reason this Appendix only refers 
to that setting. 
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Table 1: Root observation IDs of the RXTE observations used in this analysis. 



10068-08, 10122-01, 10122-02, 10122-03 

20094-01, 20095-01, 20097-01, 20415-01 

30080-01, 30081-01 

40059-01 

50136-01 

60024-01, 60025-01 

70020-01, 70020-02, 70020-03, 70021-01 
80027-01, 80114-01 
90025-01, 90426-01 



Table 2: All observations showing twin kHz QPOs with significances of more than 3a for both kHz QPOs 

(A G), and more than 3(7 for the lower kHz QPO and at least 2.5(t for the upper kHz QPO (H-K). Listed 
are the observation ID, centroid frequencies, rms amplitudes, and quality factors, as well as the centroid 
frequency difference and ratio. 



Data Obs-ID j/^ rms Q i/J rms Q Ai/ 

Set (Hz) (%) (Hz) (%) 



A 


10122-0,3-08-00 


56.1±1.3 


3.88±0.11 


0.8±0.1 


229±18 


1.41±0.24 


2.1±1.2 


173±18 


4.1±0.3 


B 


10122-01-15-01 


83.1±2.5 


2.94±0.15 


1.3±0.2 


281±32 


1.69±0.32 


1.6±0.9 


198±32 


3.4±0.4 


C 


20095-01-08-00 


97.5±2.7 


3.13±0.12 


l.OiO.l 


311±20 


1.93±0.21 


1.8±0.5 


214±20 


3.2±0.2 


D 


20097-01-32-00 


100±26 


2.74±0.40 


0.8±0.5 


419±48 


2.03±0.43 


2.6±1.6 


319±54 


4.2±1.2 


E 


10122-01-18-00 


106.5±2.5 


2.93±0.13 


1.2±0.1 


359±42 


2.28±0.25 


0.9±0.4 


252±42 


3.4±0.4 


F 


10068-08-01-00 


122.3±7.2 


2.61±0.19 


0.9±0.2 


407±33 


1.74±0.24 


1.8±0.7 


284±34 


3.3±0.3 


G 


20094-01-01-030 


127.6±3.6 


2.53±0.12 


l.OibO.l 


445±20 


1.60±0.24 


1.8±0.5 


317±20 


3.5±0.2 


H 


20094-01-01-01 


136.5±5.4 


2.65±0.21 


1.1±0.2 


404±17 


2.08±0.19 


1.9±0.5 


267±18 


3.0±0.2 


I 


20095-01-11-00 


137.4±5.4 


2.19±0.39 


1.1±0.3 


438±38 


1.47±0.31 


2. Oil. 3 


301±29 


3.2±0.3 


J 


20094-01-01-020 


157±22 


2.37±0.26 


0.6±0.2 


497±61 


1.47±0.33 


1.3±0.6 


340±64 


3.2±0.6 


K 


20094-01-01-02 


226±18 


1.24±0.28 


1.8±1.0 


505±51 


1.65±0.22 


1.5±0.5 


279±54 


2.2±0.3 



12 




Frequency (Hz) 

Fig. 1. — The power spectrum and fit function of observation 20094-01-01-01 in a power times frequency 
representation. The lower and upper kHz QPOs {Li and Lu) and the low- frequency QPO (Llf) are marked, 
as well as are the Low Frequency Noise (break component) and Very Low Frequency Noise (power law at low 
frequency). The (unmarked) component present above the low frequency QPO (here at ^22 Hz) is usually 
consistent with being a harmonic to Llf, particularly when it has high Q, cf. Fig. 8. 
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Frequency (Hz) 



Fig. 2. — Power spectra and fit functions in the power spectral density times frequency representation for 
the 11 observations showing twin kHz QPOs. The curves mark the individual Lorentzian components of the 
fit. 
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Soft Color (Crab) 

Fig. 3. — The color-color diagram of all observations of Cir X-1 that we used in 16-s time resolution (grey) 
with the data already reported by Shirey et al. (1999) plotted in a darker shade. Superimposed are the 
average colors for each of the 80 observations with kHz QPOs {white and black); see plot for definition of 
the symbols. Hard color bands are 8.5-18 keV/4.8-8.5 keV, soft color bands 3.4-4.8 keV/2-3.4 keV. 
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Fig. 4. — Hardness-intensity diagram. All symbols have the same meaning as in Fig. 3. Broad color bands 
are 6.3-13/2-6.3 keV, intensity is 2-18 keV count rate. 
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Fig. 5. — Long term light curve of Cir X-1. Each data point represents the average of one of our 497 
observations. Apart from that, all symbols have the same meaning as in Fig. 3. 
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Fig. 6. — All measured characteristic frequencies as a function of low-frequency QPO frequency vlf- Plotted 
from top to bottom are the components identified as: i„ (squares), Lg [large filled circles), the component 
{open circles) that is usually dominated by the harmonic to L^p, Llf [small filled circles), and the low 
frequency noise or break component [triangles). The 11 lower kHz QPO points that appear as one of a twin 
kHz QPO are plotted with a darker shade, and the 4 upper kHz QPOs that are only 2.5-3c7 significant with 
a lighter- centered square. Error bars are plotted but are often smaller than the symbols. 
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Fig. 7. — The fractional rms amplitude of each kHz QPO as a function of its own characteristic frequency. 
Single and twin lower and upper kHz QPOs are plotted with different symbols as indicated in the figure. 
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Fig. 8. — The ratio of the centroid frequency of the component above Llf to the centroid frequency of 
Llf vs. the coherence of the former. Clearly, when the coherence is good, this component is at twice the 
frequency oi L^p- 
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Fig. 9. — The centroid frequency of the upper kHz QPO and the LF QPO vs. that of the lower kHz QPO 
for Cir X-1 (black triangles) and other sources used by Psaltis, BeUoni & van der Khs (1999) (open circles). 
The line is a power law fit to all upper kHz QPO points. 
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Fig. 10. — Characteristic frequencies of all components vs. the characteristic frequency of L„. Previous 
observations of Z and atoll sources are plotted in black as indicated in the figure; the tracks of L^, Liowi Lh 
and Lh as previously identified in atoll sources are indicated. Colors distinguish the different components we 
identified in Cir X-1 (cf. Figs. 1 and 6): lower kHz QPO {blue), LF QPO harmonic [red), LF QPO (orange), 
low- frequency noise/break component (green). 
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Fig. 11. — The centroid frequency separation /S.v of the twin kHz QPOs vs. tlie upper kHz QPO ccntroid 
frequency. Also plotted are the straight line {dashed) and the curve for the relativistic precession model 
[drawn) best fitting the Cir X-1 data. 
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Fig. 12. — Two c'xamples of observations whose Poisson noise spectrum would have been estimated less 
accurately with other methods. 

Left: The Poisson spectrum according to the model of Zhang et al. (1995, 1996) using some previously 
suggested deadtimc values {td=S.5 /iscc, fyLB=150 /iscc) (dotted) fails to precisely describe the power at 
high frequencies. With a shift (Klein- Wolt et al. 2004) to match the > O-TbfNyquist region (dashed), power at 
frequencies up to 2 kHz might be overestimated. The Poisson spectrum computed with the deadtime values 
adopted here and no shifting (drawn) follows the data better in the range > 1600 Hz where no contribution 
from the source is expected, and is more conservative with respect to kHz features than the shifted power 
spectrum. 

Right: The strong frequency dependence of the Poisson-noisc power at very high frequencies (observed for 
high time resolution) cannot be described with the Zhang function (dotted) for the deadtime values mentioned 
above. The curve lies systematically above the data for several kHz, or partially above and partially below it 
when shifted (dashed), at frequencies where the data should purely be described with Poisson noise. Again 
the estimate with our adopted deadtime parameters (drawn) is better and appears slightly conservative with 
respect to possible high-frequency features. 
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